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1 Introduction 



The study of hypersurfaces of discontinuity in general relativity has begun 
in early twenties [1-3]. But it has been revived through new questions raised 
in cosmology and black-hole physics. Domain walls separating two coexist- 
ing different phases in inflationary scenarios [4], bubble dynamics [5], worm- 
holes [6], signature changes [7], and interior structures of black-holes [8] are 
just some of the recent applications of the thin shell formalism of general 
relativity. 

The traditional and mostly used method of handling such problems is that 
of Darmois-lsrael (DI), based on the Gauss-Kodazzi decomposition of space- 
time[9, 10]. It expresses the surface properties in terms of the jump of ex- 
trinsic curvature across the layer directly as functions of the layer's intrinsic 
coordinates. Thus the four-dimensional coordinates may be chosen freely 
and independently, adapted to the symmetry requirements, on the two sides 
of the layer. This is the very practical advantage of that method which has 
found its final formulation in the outstanding paper of Israel[9] . The geomet- 
ric conditions for the layer to be considered as a boundary of two different 
manifolds glued together at this boundary are first formulated by Darmois[3]. 
Those are minimum conditions which has been assumed by Israel. There are 
other conditions formulated by Lichnerowicz[ll, 12], which means basically 
continuous coordinates across the layer, and seems to be necessary for us- 
ing distributional tensor calculus. It is interesting to note that Sen[2], in 
this relatively unknown paper, uses the same conditions, without any fur- 
ther discusion, derives the reduced Einstein's equations for the general case, 
and solve it for the spheric symmetric 2-1-1 dimensional mass distribution. 
The O'Brien- Synge conditions [13, 14] are in most cases equivalent to Lich- 
nerowicz ones, so we are not going to consider them here[15]. Because of the 
restrictive choice of coordinates the Lichnerowicz conditions are not usually 
used. But there are cases where the distributional method, and therefore the 
Lichnerowicz conditions, are more suitable for calculation. The case of cylin- 
drically symmetric thin layer, for example, has been solved by distributional 
method using these conditions [16]. 

There have been attempts to formulate the problem of thin layer in gen- 
eral relativity using distributional methods, familiar in other part of physics. 
In fact many physical systems, classical and relativistic, undergo very rapid 
transitions of their state of motion. Think of shock waves in hydrodynam- 
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ics. Although the state of the system need not be described by discontinuous 
functions of space and time, or by functions having discontinuity in their first 
or second derivatives, but a mathematical description of the system which is 
based on distribution valued states of the system give an accurate picture of 
some important aspects of the physical problem. Usually such a description 
is more amenable to treatment than the treatment which contains a smooth 
description of the physical state. This has not only been used in classical hy- 
drodynamics, but also bas been applied in the relativistic case. Lichnerowicz 
[11, 12] has given a discussion of hydrodynamic and gravitational shock wave 
problems by using curvature tensors for space-time which contain a Dirac S 
function with support on a submanifold. Y.Choquet-Bruhat has used similar 
methods to treat high-frequency gravitational waves [17]. 

Rapid changes of physical quantities occur also in electromagnetism. One 
might have, for example, some charge distribution which is confined to a one- 
or two-dimensional region of space small compared with characteristic dis- 
tances of the problem. This distribution of charge can be replaced by a 
concentrated source, and the problem can be formulated in the sense of dis- 
tributions, and not of smooth functions. There is a natural mathematical 
framework in electromagnetism. Recall that linear operations, including dif- 
ferentiations, make sense when applied to distributions. Hence, Maxwell's 
equations, by virtue of their linearity in both fields and sources, make sense 
as equations on these distributions. This means that the machinery of distri- 
bution theory is available in electromagnetism, and guarantee that distribu- 
tional Maxwell fields with distributional charge- currents make physical sense, 
and at the same time gives a well defined and detailed sense in which a dis- 
tributional charge density must approximate our real smooth charge density 
distribution. 

In general relativity, where the field equations are non-linear, the use of 
distributional objects seem not to be trivial. Although Raju[18] claims to 

give an analytical formalism to deal with the occurrence of jump discon- 
tinuities in the metric across a hypersurface, using a non-linear theory of 
distributions [18], his method does not seem to be conclusive. The applica- 
tion he mentions is for a continuous metric where no non-linear distributional 
operations are needed. Anyhow, here the mathematical framework cannot be 
as simple as for electromagnetism. Efforts to implement distributional meth- 
ods in general relativity goes back to works of Papapetrou and Treder[20]. 
Nariai[21], aiming to understand the O'Brien-Synge junction conditions, and 
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demanding the Einstein's tensor to be free of 6 functions, uses continuous 
metrics and so brings in distributions in the formahsm, which is consequently 
used by Kumar[22]. Papapetrou and Hamoui[23, 24] then try to formulate 
a general method with application to spherical symmetric thin layer. Their 
method has then been reformulated and corrected by Evans[25]. Lichnerow- 
icz [12] gives a detailed mathematical analysis of distributions in curved 
space-time and comes to the conculsion that the classical properties of the 
covariant derivatives and all of the corresponding formulas are valid for ten- 
sor distributions. Barrabes [26] uses tensor distributions specially to include 
null hyper-surfaces in the shell dynamics. 

Taub[27] is interested in relativistic hydrodynamics and shock waves but also 
discusses the previous accounts on the concentrated 2- and 3-dimensional 
sources. Israel [2 8] and Taub[27] give a formulation for a 2-dimensional con- 
centrated mass distribution. These attempts have been criticized by Geroch 
and Traschen[28] who give an extensive and thorough analysis of concen- 
trated mass distributions in general relativity. This work is a milestone in 
all the discussions about the validity of distributional Einstein's field equa- 
tions and its applications to concentrated sources. There the authors define 
some regularity conditions for metrics, for which the distributional tensor 
calculations are allowed. Hence, for example, the line source case should be 
handled with care, although some authors have criticized it [30]. The (2+1)- 
dimensional result of this work, should not arise any problem, 

as far as the continuity of the metric is assured. We have thereafter used 
distributional tensors to solve the Einstein's equations directly, without any 
use of Gauss-Kodazzi decompositions. Although the coordinates have to be 
prepared to make them continuous at the hypersurface of discontinuity, but 
the method has been applied easily for several cases[16, 31, 32]. Naturally, 
there should be no difference in the results using either of the distributional- 
or Gauss-Kodazzi-method. But the complete equivalence has never been 
demonstrated explicitly, and the role of jump conditions by using the distri- 
butional method has never been clearly stated. 

With the results of Geroch and Traschen in mind, we show here that all 
the dynamical- and constraint-equations derived by the Dl-formalism results 
very naturally in the distributional method, without any needs to define a 
new covariant derivative. 

In section 2 we review shortly the D-I formalism and give the necessary for- 
mulae. Section III begins first with the formulation of the distributional 
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method. Covariant derivative of distributions and some useful formulae are 
given in section 3.1, and the Einstein's equations for a thin shell in section 
3.2. Section 3.3 deals with the conservation laws and the Bianchi identities. 
In this section we will see the full equivalence of the two methods. We end 
with a conclusion in section 4. 
Conventions and definitions: 

We use the signature ( — h ++), and follow the curvature conventions of 
Misner, Thorn, and Wheeler (MTW)[10]. However, our sign convention for 
extrinsic curvature is that of Israel[9], which is the opposite of MTW. The 
greek indices rum from to 3 and latin indices from 1 to 3. A semicolon 
indicates covariant derivatives with respect to either the four- metric of the 
whole space-time or to the three-metric of the layer. There will, however, 
be no confusion because the kind of indices and objects used makes the dif- 
ference transparent. The symbol denotes the covariant derivative with 
respect to either of the metrics of partial manifolds which are to be glued 
together. 

The square brackets [F] are used to indicate the jump of any quantity F at 
the layer, and bars F the arithmetic mean of it. As we are going to work 
with distributional valued tensors, there may be terms in a tensor quantity 
F proportional to some 5-f unction. These terms are indicated by F. 



2 Darmois-Isreal Formalism 

Assume two space-times M"*" and M~ with boundaries E"*" and S~. M"*" and 
may have been cut from space-times Mi and M2, respectively, but this 
is irrelevant for our task of glueing these together. Coordinates on the two 
space-time manifolds are defined independently as x+ and x'^ , and the metrics 
denoted by g''J^^(x+) and g~'^p{xt.). The induced metrics on the boundaries 
are called (y'|j({^) and g^jiC'^), where are intrinsic coordinates on E±, re- 
spectively. Bringing these 3- and 4-dimensional quantities in connection is 
trivially done with the help of tetrads defining on E[34]. 
Now, to paste the manifolds together we demand the boundaries to be iso- 
metric having the same coordinates 
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The identification 

S+ = S_ =: E 
gives us the single glued manifold M = M+ IJ M_ . 

This is the minimum requirement for glueing two manifolds together. For- 
mulated as 

[9ij] = (1) 
gives together with the continuity of the second fundamental form on E 

[K,,] = (2) 

the Darmois conditions. Both conditions should be satisfied if E is just a 
boundary surface. But in case of a thin shell we do not expect the second 
condition to be satisfied. In fact, the matter content of the shell should lead 
to a jump in the extrinsic curvature Kij. 

The condition (1) leaves the coordinates in free. If we assume the con- 
tinuity of the coordinates x± at E we then have to require 

[g,..] = 0, (3) 

which together with the corresponding equation for derivatives of the metric 

[^] = (4) 

gives the Lichnerowicz conditions. In the following we just assume the con- 
dition (1) or (3) respectively. 
On E we define a three-bein 

- A 

having the components 

The induced metric on E is given by the scalar product 

9ij = ei ■ Cj = g^uefe^j (6) 

Note that, because of the assumed isometry, this metric is the same on both 
faces E_|_ and E_. We note that the subscripts of the three-beins on E are 
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not the component indices, but as this distinction is trivial we prefer for the 
sake of simphcity not to use parentheses to distinguish them, as is usually 
done. 

We choose the parametric equation for E in the form 

Hx^ie)) = 0, (7) 

having the unit normal four- vector given by 

= a-'d^<i>, (8) 

where 



a^±J{\g-^-^^\). (9) 

Therefore 

= 0, (10) 

and 

n^n^ = e, (11) 

where e = +1 or —1 for S to be time- or space- like, respectively. We suppose 
n'^ directed from M~ to M+, i.e. in the direction of increasing a space- or 
time-like coordinate corresponding to time- or space-like S. Therefore we 
have to take the positive (negative) sign in (9) for time- (space-)like E. This 
choice gives us the useful relation 

lal 

sign a = — = e. 
a 

The choise of Lichnerowicz condition (3) makes it possible to have a unique 
normal vector for each case. As wc want to concentrate on the formulation 
of the distributional method and avoid any undue complications, we leave 
aside the case of null hypersurfaces. 

Now, in general the metrics in M+ and M~ need not to be continuous at E, 
but they could be. However, the normal extrinsic curvature (second funda- 
mental form) is not continuous for a thin shell. It is defined by 

= -n^e^.V^.e';. = (12) 
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Now, we have all the prerequisites to write the Einstein's equation for the 
hypersurface. These are 10 equations which will be written in components 
normal and tangent to the hypersurface. The first and second contracted 
Gauss-Kodazzi equations are [9, 10] 

G^.rfn'' = ^{K^ - K^jK'^ - e 'R) (13) 

G^.e'ln'' = K\.^ - K„ (14) 

where % and \} are the Ricci scalar and Einstein tensor of the three metric 
Qij, respectively. Now, to discover the effect the energy- momentum tensor 
Sij of E on the space-time geometry, we perform a "pill-box" integration of 
Einstein's equations across E: 

S^^ = lim (^T^^ ~ ^"""k) K iH?o i_s ^^^-^ 

where n is the proper distance through E in the direction of the normal n^. 
Sfji_y is the associated 4-tensor of energy momentum of the shell. The equa- 
tions (13, 14) have the physical meaning that no moment associated with the 
surface layer flows out of E. Therefore S^j, vanishes off the hypersurface E, 
which is expressed as 

S^^rf = (16) 

The energy momentum 4- and 3-tensors are related as 

Si"^ = ete'^S'^ (17) 

The covariant derivative of such a tensor relative to the corresponding con- 
nections is given by [9] 

V^S'^^ = ef - eS'^K%n^, (18) 

which leads to the following useful relation 

ey.Sl = Si^ (19) 

Similarly we can associate to the 3-tensor Kij defined on E, the corresponding 
4-dimensional tensor: 

Kt^y ^ K'^efe'^j, (20) 

satisfying 

K^'^n^ = 0. (21) 

The remaining components of the Einstein's equations lead to the following 
non- vanishing result 
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lim / Guueie", dn = 

emj]-g,j[K])^KS,j (22) 

This distributional equivalent of Einstein's equations is called Lanczos equa- 
tion, which partly determines the dynamic of the thin shell. The other dy- 
namical equations come from the defining equation of matter contents of 
the shell. Now, the two Gauss-Kodazzi equations act as constraints. The 
first one (13) is the so-called "Hamiltonian "- and the second one (14) the 
"ADM"- constraint. Note however that these equations are vahd in M+ and 
M~ on taking the limits as one approaches the layer E. Therefore we are 
actually faced with 8 equations, the sum and difference of which give us the 
junction conditions. The Hamiltonian constraint along with the Einstein's 
or Lanczos equations then give the evolution identity: 

S'^Kij^-iT^^n^'n'' -A/k] (23) 

and 

+ (KijTt' - K^) = 2eK{T^,nf^n- - A/k) + ^-^{S'^ - ^ ). (24) 
The ADM constraint gives the conservation identity 

S% = -efT^.n'^eJ] (25) 

and 

X'^.,-Z,=«:(7>^ (26) 

Not all of these jump conditions are independent. Usually one takes the evo- 
lution identity (23) and the conservation identity (25) as the proper junction 
conditions, which in addition to the Lanczos equation should be satisfied [35]. 

3 Distributional Method 

Here we intend to give a formulation of the Einstein's equations for the 
case where there exists a hypersurface of concentrated source immersed in 
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an otherwise arbitrary space-time, not necessarily vacuum. We assume the 
metric to be continuous at the hypersurface: 

M = (27) 

Otherwise we would have to consider non-linear operations of distributions 
such as 56 or 55. The disadvantage of having a continuous metric across the 
shell pays off by the simplicity of the method to calculate specific solutions[16, 
31, 32]. 

Write the metric in the following form 

g,.^glAm)+9-,A-^{x)) (28) 
where 6 is the step function and 

g%W=o)=G = g~fxvW^)=G (29) 

This condition guarantee the smoothness of the metric on the hypersurface. 
Should this not be the case we try a coordinate transformation x — x{x') 
having a jump in the first derivative: 

^ = a%^e{^{x)) + a-^e{-^{x)) (30) 
The condition for the new metric to be continuous comes out to be 

We assume from now on that the metric is smooth everywhere, at the 
hypersurface and on both sides of it. 

Although the metric is continuous on S, its derivatives, and so the corre- 
sponding connections, are discontinuous. Nevertheless the connection corre- 
sponding to the metric g^j^v can be written in the following compact form: 

= ^($(a;))r+V + ^(-*(^))r-V' (32) 

where are the ordinary connections on M^. The above connection has 

jump discontinuities on E. 
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To write the field equations for tlie liypersurface we need tlie formulation of 
the energy-momentum tensor of the shell. Generally it can be written in the 
form 

= CS^,5{^{x)) (33) 

where C is a constant to be calculated. We integrate the above equation in 
the direction of the normal to the hypersurface 

j %,dn^CS^, j 5{^{x))dn^CS^,\^\ (34) 
Therefore, using the definition (15), we obtain 

-Hfl (3«) 

and 

%u^S,,\'^\8{^{x)) (36) 

Note that in the hterature one usually takes C — 1, which is correct just for 
special cases. But in general the factor C is necessary (see also [33]). Now, 
the derivative of $ in the normal direction can be written in terms of unit 
normal vector n'*: 

C = I — I = = leal = lal, (37) 

an 

where we have used (8-11). Therefore (33) will be written in the form 

= CS^,5{^x)) = \a\S^,5{^{x)) (38) 



3.1 Covariant Derivative of Distributional valued Ten- 
sors 

There is no need to change of the ordinary concept of covariant derivative, 
as it has been carefully shown by Lichnerowicz [12]. In fact, all the known 
properties of covariant derivative and the corresponding formulae in a pseudo- 
Riemannian manifold are valid for tensor distributions. But for the sake of 
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convenience of calculation we refer to some useful formulae. Consider first 
an arbitrary vector A'^ defined as 

A^" ^e{^)A+ + e{-^)A-, (39) 

where A"^ has the support on M^. It is therefore useful to define the operator 

V = ^($)V+ + ^(-$)V- (40) 

We can know write the covariant derivative of a distributional valued vector 
A^^ in terms of the covariant derivatives of its defining parts in M^. The 
following relation is easily obtained: 

Ai^^^ = V.A'^ + [A^^] 9,$ 5($) (41) 

This relation can be generalized easily for a distributional tensor of any rank. 
The covariant derivative of the tensor 

T^p) = e{^)T+^f^ + e{-^)T-^P\ (42) 

where (p) stands for any number of indices, is calculated to be 

T^p).^ = V^T^p'^ + [T^"^] 9,$ 5($) (43) 

In the case a tensor has the support on E its covariant derivative is in the 
usual form. Take the tensor T^''^ from (33). Its covariant derivative can be 
written 

{f>"'),p = (C5^");p 6{<^) + CS"" (5($)),p 

= {CS^"'),p6{^) + CS^"' dp<i>6'{^) (44) 
We will need this relation later to discuss the conservation laws. 

3.2 The Field Equations 

The Einstein's field equations are valid on both sides of the hypersurface as 
usual. So we concentrate our procedure on E, where we expect the curvature 
and Einstein tensor to be proportional to 5. That means in calculating the 
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connection coefficients and the components of the Ricci tensor we can ignore 
terms not proportional to 6. Hence, e.g., the terms in the Ricci tensor 

— ^fip,u ~ ^'^lu,p ^i^p^au ~ r^^r^^ (45) 

proportional to F's can be ignored. The only relevant terms are 
Now, 

r^P = 4^,.. (47) 

where g is the determinant of the metric. The 5 distribution can only occur 
in the second derivatives of the metric. Therefore 

K,^ = (48) 

Similarly, for the second term in the Ricci tensor we have 

1 ^ ^ 

tii',p ~ 2^'' {9(Tfx,vp + 9(yv,fxp ~ 9ixv,ap) (49) 

Having the metric in the form (28) we obtain 

gaP,,.^[gapA{dr^^)^{^{^)) (50) 

and 

l,u = [g^dMi^ix)) (51) 
As the result we obtain for terms in the Ricci tensor proportional to 5 

R,u - [^^gAd.^-riM + lg^^A-l^^A)^!^^^ 



[^J9.$-[rM9,$U($(x)) (52) 



This enable us to write the Einstein's equations for the layer: 

Gp,!/ = kT^^ (53) 
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Defining 




- i^Uj (54) 

we obtain, using (38) and (52) for tlie energy momentum tensor tlie field 
equations in the 4-dimensional form 

Qfiu ~ i^9iJii^Q = ^i^S^y, (55) 

where Q — Q^vg^", and we have used the relation e = Note that Q^i, is 
a tensor with support on E. This equation, for the time-like case, has been 
first derived, without using the hitherto unknown distributional calculus, by 
Sen[2]. We would like, therefore, to coin it by Sen equation. The three 
dimensional form of the Sen equation is readily obtained by decomposing it 
to tangential- and normal-components to E. Multiplying (54) with we 
obtain 

S^uu" = 0, (56) 

which is the same relation as (16). This tell us immediately that the compo- 
nents corresponding to Sn^n'^n^ and S^^n^e'^ identically vanishes. To obtain 
the proper 3-dimensional components we notice first that 

Q^J = Q,ue^e^ = -KM,e>te) = [K,,]. (57) 
Threfore, we obtain from the Sen equation 

Qij = eK (^Sij - ^9ijS^ , (58) 

which is equivalent to the Lanczos equation (22). 

We have therefore seen that the exphcit method of writing the Einstein's field 
equations for a regular metric which is continuous without having continuous 
derivatives leads to the equation (55) and is equivalent to the Dl-formalism 
based on the Gauss-Codazzi formalism. In practice, one begins with known 
solutions of the Einstein's equations in M"^, and after making sure the conti- 
nuity of the metric on E, tries to solve the equations (55). In the following we 
will show that the jump conditions of Dl-formalism follows from the Bianchi 
identities corresponding to the metric (28), and are therefore implicit in the 
equations we have used. 
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3.3 Conservation Laws 



We have now all the prerequisites to evaluate the Bianchi identities and the 
conservation of energy momentum tensor of our pasted space-time. The 
energy momentum of the whole space-time, including the cosmological terms 
A± is 

rpixu ^ ft.u ^ _ ^+ 1^ g+f'^)e{^) + (T-/^^ - K-/k g-^"')e{-<^>) , (59) 

where T^'^ is defined in (38). Having in mind that the covariant divergences 
of T^^*^ and g'^^''' with respect to the corresponding connections vanishes, we 
obtain 

T^- ^ = (f '^-).^ + [T^^- - A/n] 5($), (60) 

where we have used the relation (43). Now inserting for T^'^ from (38) and 
using (44) we obtain an equation having terms proportional to 5($) and 
5'($). Each term vanishes independently. The term proportional to 5'($) 
gives 

{S'"'d,^){n''da^) ^0, (61) 

which leads to the (56) and ensures the orthogonality of the energy-momentum 
tensor of E to the hypersurface normal n'^. We use in the following this re- 
lation to simplify the remaining calculations. The term proportional to 5 
is 

[CS^""),^ + CS^'Prip + CSPT^p^ + +[T^'' - A/k g^^'jd^^) 5($) = 0. (62) 

Wc have left 5($) as proportionality factor to stress its influence specially 
on terms containing F's and g^'^^s. Note that the third term containing F's 
contains terms like 9 ■ 5, i. e. product of distributions. This is in analogy 
to the elementary problem of evaluating the electrostatic force on a sheet of 
charge [5]. There the linearity of electrostatic equation resolve the ambiguity. 
But how about our case where the Einstein equations are non-linear? We 
have already shown that in the case of thin shells the only terms contributing 
to the Einstein tensor are the derivetives of the connection, or the second 
derivative of the metric, which appear linearly. It is then easily seen that 
is case of a layer the Einstein's equations leads to a Poisson-like equation 
corresponding to the Sen equation (55) for concentrated distribution, where 



15 



the second derivative can be replaced by (50) and (51). Therefore, in analogy 
to the electromagnetic case we can use the linearity of (55) to show that 

r^.<^(^) = ^(r+V + r-V)<^(^) (63) 

Using this result and multiplying the equation (62) with n^, we obtain the 
component in the normal direction. 

{CSn,^n, + C5^^r%n^ = e[T^^ - A/k g'''']n^d,^. (64) 

Using the relations (63) and the definition of the extrinsic curvature (12), we 
obtain the final result 

S'^Kij = efT'^'^n^n^ - A/k] (65) 

This is the evolution identity (23) derived as one of the jump conditions 
in the Darmois- Israel method. Here it is just a consequence of the Bianchi 
identities. To obtain the remaining three equation we multiply (64) with e^. 
Using (19) we obtain 

S)^, = -e [T,,n^e\] , (66) 

which is the conservation identity (25). It gives the conservation law for the 
energy momentum tensor of the layer. We therefore see that our explicit 
distributional method of solving the Einstein's equations gives all the dy- 
namical and constraint equations of Darmois-Israel method, and is therefore 
equivalent to it. 



4 Conclusion 

We have seen that, based on the Lichnerowicz condition (3), a distributional 

method can be formulated to solve the Einstein's field equations for a thin 
shell, which is equivalent to the Darmois-Israel formalism, and gives all the 
necessary equations and jump conditions formulated there. In fact, it has 
been shown that the jump conditions are consequence of the Bianchi identi- 
ties, and therefore implicit in the formalism, once the Lichnerowicz condition 
is satisfied. This makes the distributional formalism easy to apply, specially 
when explicit solutions are to be found, and it pays off the disadvantage of 
the continuity of the coordinates across the shell. 
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